riep-ph/01 11209 



FERMILAB-Pub-01 /348-T 



Two-loop effective potential for a general renormalizable theory 
and softly broken supersymmetry 



Stephen P. Martin 

Department of Physics, Northern Illinois University, DeKalb IL 60115 and 
Fermi National Accelerator Laboratory, P.O. Box 500, Batavia IL 60510 



Abstract 

I compute the two-loop effective potential in the Landau gauge for a general renormal- 
izable field theory in four dimensions. Results are presented for the MS renormalization 
scheme based on dimensional regularization, and for the DR and DR schemes based 
on regularization by dimensional reduction. The last of these is appropriate for mod- 
els with softly broken supersymmetry, such as the Minimal Supersymmetric Standard 
Model. I find the parameter redefinition which relates the DR and DR schemes at 
two-loop order. I also discuss the renormalization group invariance of the two-loop 
effective potential, and compute the anomalous dimensions for scalars and the beta 
function for the vacuum energy at two-loop order in softly broken supersymmetry. 
Several illustrative examples and consistency checks are included. 
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1 Introduction 

The Fermilab Tevatron collider and the CERN LHC collider hold the promise of exposing 
the nature of spontaneous electroweak symmetry breaking. In the Standard Model, this 
mechanism relies on a non-zero vacuum expectation value (VEV) for a fundamental Higgs 
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scalar field. There are good theoretical and experimental reasons to suspect that this picture 
is correct, but incomplete, and must be embedded in a larger theory such as supersymmetry 
0, When new experimental discoveries are made, the tasks of telling the difference 
between different candidate models of electroweak symmetry breaking and constraining the 
underlying parameters of the successful theory will require high-precision calculational tools 
at the two-loop level or better. 

The effective potential allows the calculation of the VEVs in the true vacuum state 

of a theory with spontaneous symmetry breaking. In this formalism, the scalar fields of the 
theory are each separated into a constant classical background <fi plus quantum fluctuations. 
The effective potential V(4>) is equal to the tree-level potential in the classical background, 
plus the sum of one-particle-irreducible connected vacuum graphs. These are calculated 
using the Feynman rules with ^-dependent masses and couplings. Thus one may write 

v = v m + A_ vm + _L_ v v + ..., (L1) 

where represents the n-loop correction.^] In this paper, I will be concerned with the ef- 
fective potential in Landau gauge. Although the effective potential itself is gauge-dependent, 
physical properties following from it, such as its value at stationary points, and the ques- 
tion of whether or not spontaneous symmetry breaking occurs, are gauge invariant . The 
one-loop contribution is well-known for a general field theory, and is reviewed in section 
|3|. In ref. |7|], Ford, Jack and Jones have calculated in the special case of the Standard 
Model using dimensional regularization (DREG) with minimal subtraction or modified min- 
imal subtraction (MS). Their calculations can be generalized to obtain the corresponding 
result for any renormalizable field theory, as I will do explicitly in section ^. 

However, it is well-known that the DREG regularization method is not convenient for 
theories based on supersymmetry. This is because in DREG, the vector field only has 
4 — 2e components, introducing a spurious non-supersymmetric mismatch with the number 
of degrees of freedom of the gaugino. Therefore, in DREG the relationships between couplings 
which should hold in a softly broken supersymmetric theory are violated even at one-loop 
order. Instead, one can use the dimensional reduction (DRED) method |§, in which loop 
integrals are still regularized by taking momenta in 4 — 2e dimensions, but all 4 components 
of each vector field are kept. The extra 2e components of the gauge field in DRED transform 
like scalars in the adjoint representation of the gauge group, and are known as epsilon 

^To save ink, a factor of l/(167r 2 )™ is always factored out of the n-loop contribution to the loop expansion 
of the effective potential, /3-functions, and anomalous dimensions in this paper. 
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scalars. The renormalization scheme based on DRED with modified minimal subtraction is 
known as DR. It has the virtue of maintaining manifest supersymmetry in theories where 
supersymmetry is not explicitly broken. 

Realistic models of the physics at the TeV scale do involve explicit soft violations of 



supersymmetry, however. In such models, the DR renormalized dimensionless couplings of 
the theory obey the relations prescribed by unbroken supersymmetry. However, the epsilon 
scalars in general do not have the same masses or dimensionful couplings as do the ordinary 
4 — 2e vector field. In fact, computation of the renormalization group (RG) equations shows 
that the running squared masses of the epsilon scalars cannot be consistently set equal 
to those of the corresponding vector gauge bosons f9j. This makes the DR scheme also 



inconvenient, since the epsilon-scalar masses are unphysical. A better scheme is the DR 



scheme [10[ , which differs from DR by a parameter redefinition. The DR scheme offers the 
advantages that the epsilon-scalar masses completely decouple from all RG equations, and 
also from the equations that relate running renormalized parameters to pole masses and 
other physical observables. 

In this paper, I will present results for the two-loop effective potential in the Landau 
gauge and in each of the MS, DR, and DR' renormalization schemes. For models with exact 
supersymmetry, the last two schemes are the same, while for models with softly broken 



supersymmetry the DR scheme is by far the most convenient. 

The topologies of the one-particle-irreducible connected vacuum graphs at one- and two- 
loop orders are shown in Figure [I| Because the one-loop graph topology does not involve 
interaction vertices, clearly depends only on the field-dependent squared masses m 2 , 
where the index n runs over all of the real scalars, two-component fermions, and vector 
degrees of freedom in the theory. Note that any complex scalar can be written in terms of 
two real scalars, while four- component Dirac and Majorana fermions can always be written 
in terms of two-component left-handed Weyl fermions, in a way throughly familiar to dis- 
ciples of supersymmetry (see refs. [|IJ |2|] for a discussion). In any dimensional-continuation 
regularization scheme, quadratic divergences are automatically discarded, and one finds for 
the renormalized effective potential at one-loop order: 

V® = - £(-l) 2 -(2 Sn + l)(m 2 ) 2 (hrm 2 - k n ). (1.2) 

^ n 

Here I have adopted the notation 

E(x) = ln(x/Q 2 ), (1.3) 
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Figure 1: Topologies of one-particle-irreducible connected vacuum Feynman diagrams for 
the one-loop and two-loop contributions to the effective potential. 

where Q is the renormalization scale, and s n = 0,1/2,1 for real scalars, two-component 



where g npqr and g npq are field-dependent four- and three-particle couplings, and f n p{x,y) 
and f npg (x,y, z) are Q-dependent functions obtained by performing the appropriate two- 
loop integrations. So the task is to identify these objects for each combination of particle 
types that can contribute. 

The rest of this paper is organized as follows. Section 2 describes the field- dependent 
masses and couplings, lists the relevant Feynman diagrams, and presents necessary conven- 
tions. Section 3 reviews the one-loop effective potential, distinguishing between the MS, 
DR, and DR' schemes. Sections 4-6 present the results for the two-loop effective potential 
contribution in each of those schemes. Section 6 also explicitly gives the redefinitions neces- 
sary to go from DR to DR'. Section 7 discusses the RG invariance of the effective potential 
in the DR' scheme, and derives some necessary results for the scalar anomalous dimension 
and vacuum energy beta function in softly broken supersymmetry. Section 8 contains some 
illustrative examples and consistency checks. 




2 Conventions and setup 



*The contribution of epsilon scalars is discussed in section 
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2.1 Field-dependent masses and couplings 

Let us write the quantum fields of a general renormalizable field theory as a set of real 
scalars R' t , two-component Weyl fermions ip'j, and vector fields A'£. Scalar flavor indices are 
i, j, k, . . .; fermion flavor indices are /, J,K,.. .; and a,b,c, . . . run over the adjoint represen- 
tation of the gauge group. Space-time vector indices are written as Greek letters p,,u,p,.... 
I use a metric with signature ( — h++), and the notations for fermions follow [p], [|. The 
primes are used to indicate that these fields are not squared-mass eigenstates. The kinetic 
part of the lagrangian includes 

-£= \n%R\R) + \{m IJ ^'j + c.c.) + l -ml b ^^ b . (2.1) 

The symmetric fermion mass matrix m IJ yields a fermion squared-mass matrix 

m ]j = m* IK m KJ . (2.2) 

Then m|- and m 2 ab are real symmetric matrices, and m u is a Hermitian matrix, and in 
general they all depend on the classical background scalar fields. In order to calculate the 
effective potential, the first step is to rotate to squared-mass eigenstate bases Ri, ipi, A^. 
This can be done by using orthogonal matrices , for the scalar and vector degrees 
of freedom, and a unitary matrix for the fermion degrees of freedom. So, the rotations 

R\ = NjpRj, (2.3) 
V4 = MPVj, (2.4) 
K = NiP^, (2-5) 



are chosen such that: 



N^mlNf = 5 l3 ml (2.6) 
nTMklNTJ* = Sum), (2.7) 
N^WcXT = (2.8) 

Here mf , mf, and m\ are respectively the scalar, fermion, and vector squared-mass eigen- 
values which will appear in the effective potential. It should be noted that in general 
diagonalizes the fermion squared-mass matrix mjj, but need not diagonalize the fermion 
mass matrix m IJ . All that is required is that 

M IJ = N$*m KL N!£* (2.9) 
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has a block diagonal form, with non-zero entries only between states with the same squared- 
mass eigenvalue. Indeed, it is quite often not particularly desirable for N^ F ' to diagonalize 
the fermion mass matrix, for example in the case of charged Dirac fermions, with doubly- 
degenerate eigenvalues for mf , where M IJ is best left off-diagonal in 2 x 2 blocks. The matrix 
M IJ and its complex conjugate MJj will appear as mass insertions. In practical applications, 
the diagonalizations just described are easily performed numerically using a computer, and 
under favorable circumstances (such as those studied in section ^|) they can be done analyt- 
ically. In either case, the problem amounts to finding the orthonormal eigenvectors of mf -, 
mjj, and m 2 ab . 

Now the interaction terms in a general renormalizable theory can be written in terms of 
the squared-mass eigenstate fields as 

£s = —-X^RiRjRk — -^^ kl RiRjRkRu (2.10) 

C SF = -^y IJk ipiipjR k + c.c, (2.11) 

£ S v = —g^A^Ri - ^g abij A^A^ RiRj - g^A^R^Rj, (2.12) 
C FV = gfA^a^j, (2.13) 
Asaugo = g abc A a ^A b u d>*A uc - -cfteg^A^A^A^A* + g abc A a ^ b d^uj c , (2.14) 

where u a and uJ a are massless (in Landau gauge) ghost fields. This defines the field- dependent 
couplings to be used in the two-loop effective potential calculation. The scalar interaction 
couplings X 1 ^ and \^ hl are each completely symmetric under interchange of indices, and real. 
The Yukawa couplings y IJk are symmetric under interchange of the fermion flavor indices 
J, J. The remaining couplings all have their origins in gauge interactions. The vector-scalar- 
scalar coupling g at i is antisymmetric under interchange of The pure gauge interaction 

g abc ig 

completely antisymmetric; it is determined by the original gauge coupling g, the 
antisymmetric structure constants f abc of the gauge group, and N^ v \ according to 

g ^ = gr f 'NSpN!pN^. (2.15) 



Similarly, if the fermions ip'i transform under the gauge group with representation matrices 
(T a )/, then the vector-fermion- fermion couplings are 

9f-9(T b ^Nf>N^N^. (2.16) 

Note that even the dimensionless couplings generically depend on the classical scalar back- 
ground fields 0, through their dependence on the rotation matrices N^ s \ N^ F \ and N^ v \ 
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2.2 The Feynman diagrams 

The two-loop effective potential is to be evaluated by computing the one-particle-irreducible 
connected vacuum Feynman diagrams appearing in figure |2|. The masses and couplings of 
fields appearing in these diagrams are as indicated on the right sides of eqs. (|2.6| )-( ]04l ). 
Dashed lines denote scalar propagators. Solid lines denote fermion propagators with helicity 
along the direction of the arrow, and large dots between opposing arrows denote insertions 
of the fermion mass matrix M IJ or its complex conjugate Mjj, depending on whether the 
arrows are incoming or outgoing. Vector propagators are indicated by wavy lines, and ghost 
propagators by dotted lines. Each graph is also labelled by the type of propagators it 
contains, with S,F,V,g standing respectively for scalar, fermion, vector and ghost. Also, 
the presence of mass insertions in fermion lines is indicated by the overlines in the labels 



FFS and FFV. The results for these Feynman diagrams (plus counterterms) are reported 
in sections [|, ||, and |6]. 

2.3 Two-loop integral functions needed for vacuum graphs 

All of the effective potential two-loop integrals can be expressed in terms of linear combi- 
nations of functions introduced and studied by Ford, Jack and Jones in 0. I will follow a 
notation similar but not identical to theirs: the functions I(x,y,z), J(x,y), and J(x) used 
here are equal to the e-independent parts of the functions I(x,y,z), J(x,y), and J(x) used 
in ref. [[?]], up to obvious factors of 1/167T 2 . Explicitly, I choose to express results in terms of: 

J{x) = x(Ex-l), (2.17) 
J(x,y) = xy(hix- l)(my-l), (2.18) 

I(x,y,z) = -{x — y — z)\ny\iaz + -{y — x — z)lnxlnz + ~{z — x — y)lurlm/ 

5 i 

+2x\nx + 2y\mj + 2zlnz — -{x + y + z) — ^(x, y, z) . (2.19) 

Here £(x,y,z) was originally found in terms of Lobachevskiy's function or related integrals 
in ref. [0 using methods developed in |TT], but it can also be expressed [0, [TJ|, |TJ] in 
terms of dilogarithms according to (for x,y < z): 

£(x,y,z) = R{2ln[(z + x -y - R)/2z]ln[(z + y -x - R)/2z] - \n(x / z)ln(y / z) 

-2Li 2 [(z + x-y- R)/2z\ - 2U 2 [(z + y-x- R)/2z] + vr 2 /3} (2.20) 

with 

R = [ x 2 + y 2 + z 2 - 2xy - 2xz - 2yz} 1/2 . (2.21) 
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Figure 2: The one-particle- irreducible connected Feynman diagrams contributing to the two- 
loop effective potential. Dashed lines denote real scalars, solid lines denote Weyl fermions 
carrying helicity along the arrow direction, wavy lines are for vector bosons, and dotted lines 
are for ghosts. The large dots between opposing arrows on the fermion lines in the FFS and 
FFV diagrams denote mass insertions. The FFS diagram is accompanied by its complex 
conjugate (the same diagram with all arrows reversed). 



9 



The dilogarithm function is denned in the standard |L6[ way for any complex z: 

ln(l - t) 



Li 2 z = - 



-dt. 



(2.22) 



10 t 

To resolve branch cut ambiguities which could arise, consistently choose Arg(_R) = or n/2 
along with 

- 7T < Im[ln(z)] < 7T (2.23) 

for all logarithms of negative or complex z, including the logarithm appearing in the integral 
definition of the dilogarithm. So, for example, when x is real and greater than 1, 

Im[Li 2 (a;)] = -i7rln(;r), (2.24) 
lm[U 2 (x±i8)} = ±ivrln(x), (2.25) 

for S real and infinitesimal, while Im[Li 2 (x)] = if x is real and less than 1. The functions 
£(x, y, z) and therefore I(x, y, z) are invariant under interchange of any two of x, y, z. 

It is useful to have expressions for these functions in the special cases of vanishing argu- 
ments. In addition to the trivial identities J(0) = 0, J(x, 0) = J(0, x) = and 1(0, 0, 0) = 0, 



one finds [ D| : 

1(0, x,y) 



(x-y) 



— 1 — 7T 2 

^h(y/x) - h\(x/y)hx(x -y) + -(ma;) 2 - — 

2 o 



-— (x + y) + 2xlnx + 2ylny — xlnxln?/, 



2J(x) — 2x J(x, x) 

x 



1 — 2 — 5 7T 

--xllna;) + 2xmx x x 

2 v ' 2 6 



x(lna;) 2 + 4xlnx — 5x, 

2 



1(0, x, x) 
1(0,0, x) 

It is also sometimes useful to expand these functions for infinitesimal arguments: 

1(5, x, y) = 1(0, x, y) + 6{-(x + y)I(0, x, y) - 2J(x, y) + 3xJ(x) + 3yJ(y) - y.J(x) 



(2.26) 
(2.27) 
(2.28) 



1(8, x, x) 



I(5 1 ,5 2 ,x) 



-xJ(y) -(x + yf + (x- y)[J(y) - J(x)]\n6\/(x - y) 2 + 

2J(x) — 2x J(x, x) 

x 

+5 1 4 + ^ J(x, x) + h(x) - [1 + ^]E* j + • • • , 

1(0, 0, x) + — \-x - 1(0, 0, x) + 3 J(x) - J(x)ln5i 
x 1 

X 

+4J(x 



(2.29) 



(2.30) 



x 



1(0, 0, x) + 3J(x) - J(x)E5 2 l + J_i{_2/(0, 0, x) 

J x 2 L 



2x — [x + J(x)}(ln5i + ln<5 2 ) + xln5iln5 2 | + . . . , 
where the ellipses stand for terms with more than one power of S or either #i or <5 2 . 



(2.31) 
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2.4 Conventions for softly-broken supersymmetric models 

One of the main applications of the results of this paper is to models with softly broken super- 
symmetry, such as the Minimal Supersymmetric Standard Model (MSSM). Therefore I now 
list the relevant conventions to be used here for such models. In general, the superpotential 
is given in terms of the chiral superfields by 

W = -Y^Qi^Qk + ^//'•'<IvI< ; , (2.32) 

and the soft supersymmetry-breaking part of the Lagrangian is 

- £ soft = {^ k fafafa + h^fafa + cV* + \MX a X a + cc.) + (m 2 )^*^ + A, (2.33) 

where the fa are the complex scalar field components of the <&j, and the A a are the two- 
component gaugino fermions with mass M. The parameter c l can only appear if there is a 
gauge-singlet chiral superfield in the theory. Note the presence of a vacuum energy term A. 
This is required in order for the full effective potential to be RG invariant [I7|-|[2U||. The two- 



loop beta function for A is obtained in section [I], and the beta functions for each of the other 



couplings at two-loop order are given in [ 21] , [22] , |9], |10[, Flipping the heights on all indices of 
a coupling implies complex conjugation, so = {Y^ k )* , ^ = (/i*- 2 )*, = (a w ' fc )*, etc. 
The representation matrices for the chiral superfields are denoted by (T a )l. They satisfy 

[T a ,T b }=if abc T c , (2.34) 

where f abc are the totally antisymmetric structure constants of the gauge group G. Then 

(T a T a )l = C^5{, (2-35) 
Tr [T a T b ] = S{R)5 ab , (2.36) 

jacdjbcd = fj^ab ( 237 ) 

define the quadratic Casimir invariant C(i) for each representation, the total Dynkin index 
S(R) summed over all representations, and the Casimir invariant of the adjoint representa- 
tion. The dimension of the adjoint representation is 

d G = Tr[C(i)]/S(R). (2.38) 

I use a normalization such that each fundamental representation of SU(N) has C(i) = 
(N 2 - 1)/2N and contributes 1/2 to S(R). 
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3 One-loop effective potential in the MS, DR, and DR' 
schemes 



In this section, I review the known answers for the one-loop effective potential. This will 
allow us to carefully distinguish the results in the MS, DR, and DR' schemes. 
In the MS scheme and Landau gauge, one has 

^>=y« + V<W ) (3-D 

where the different contributions arise from scalars, fermions, and vectors going around the 
loop in figure [l]: 

VP = ^K 2 f(K 2 -3/2), (3.2) 

i 

VP = -VK) 2 «-3/2), (3.3) 
z / 

VP = j£(^) 2 (Em 2 -5/6). (3.4) 

a 

The appearance of 5/6 rather than 3/2 in Vy is due to the fact that there are only 4 — 2e, 
rather than 4, vector degrees of freedom in MS. 

In the DR scheme, one must include also the effects of the epsilon scalars. Now, it is 
tempting to assume that the epsilon scalars have the same field-dependent mass as their 
4 — 2e vector counterparts. However, as pointed out in ref. f9|, this is actually inconsistent 
except in models with exact supersymmetry, unless one sticks to only one fixed value of 
the renormalization scale Q, because the epsilon-scalar squared mass has a beta function 
which is not homogeneous. Therefore, in general one must allow the epsilon scalars to have 
squared-mass eigenvalues rh 2 a which are distinct from the m 2 for the ordinary vectors. To 
be specific, consider the explicit form of the field-dependent squared-mass matrix for the 
ordinary 4 — 2e vector fields: 

m 2 ab = g a g b( j ) * l {T a ,T b y i ^. (3.5) 
This has eigenvalues m 2 . For the epsilon-scalar squared-mass matrix, one has instead: 

™lb = m lb + $abm 2 , (3.6) 



where mj is an "evanescent" parameter. This matrix requires an orthogonal diagonal- 



ization matrix which differs from N^: 

N^KXi = S*K ( 3 - 7 ) 
12 



Unless supersymmetry is not explicitly broken, the eigenvalues m 2 will in general differ 

from m 2 , and the corresponding couplings of the squared-mass eigenstate epsilon scalars are 

different from the couplings of squared-mass eigenstate vectors, because differs from 
N (v)_ 



In the DR scheme, with epsilon scalars included, one now finds 

V m = V P + V P + V v ] + V e (1) (3-8) 
where vft } £1X6 clS before, and 

K W = -^E(-a) 2 - (3.9) 

However, m 2 is an additional parameter with no physically observable counterpart, and 
so its appearance in the effective potential is quite inconvenient. The functional form of the 
effective potential is also not directly physically observable, so there is no contradiction; m 2 
must cancel only from observable quantities. However, clearly one would like to avoid having 
to include a distinct epsilon-scalar mass in calculations in the first place. This problem was 
solved in the context of softly broken supersymmetric models in ref. [0] with the introduction 



of the DR scheme. The point is that one can remove the dependence of the full one-loop 
effective potential on m\ by redefining the ordinary scalar squared masses and the vacuum 
energy term appearing in the tree- level part eq. (|2.33|) : 

Kr'H = ( m m)i-J^[ 5 i 2 9 2 C(i)mi\, (3.10) 



1 d G (m 2 ) 2 

A dr' = A DR~I^ 2 ' ^ ' 



(These hold for a simple or U (1) gauge group. If there is more than one simple or U(l) gauge 
group, then the correction terms should be summed over subgroups, with a distinct m 2 for 



each subgroup.) The result is the DR scheme, and the effective potential in this scheme is 



the one usually quoted in the literature (and often slightly incorrectly referred to as the DR 
one): 



v m! = E(- 1 ) 2s "( 2 ^ + 1 )MO = STr[/ i K)j, (3.12) 

n 

where 

2 

h(x) = ?j [E(ar) - 3/2] , (3.13) 
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and n runs over all real scalar, Weyl fermion, and vector degrees of freedom. The scalar 
squared masses occurring in eq. ( |3.12| ) are the ones following from the redefinition in eq. ( ft.lOQ , 



and the vector squared masses are the eigenvalues of eq. (3J5). The DR effective potential 
is both manifestly supersymmetric when the soft terms vanish, and independent of the un- 
physical evanescent parameter m 2 when the soft terms do not vanish. It is not hard to see 
that m 2 is simultaneously banished from the equations which relate the physical pole masses 
to the Q-dependent running masses in the theory, so it has been successfully decoupled from 



all practical calculations. It would be quite clumsy to use the original DR scheme in stud- 
ies of realistic models like the MSSM, since in RG running and evaluation of the effective 
potential one would have to keep extra contributions from epsilon-scalar masses in order to 
avoid inconsistencies. Therefore the DR' scheme is the preferred one. 



After making this painful distinction, it must be admitted that the DR final result for 
the effective potential has exactly the same form that one would have obtained if one had 
naively set m 2 equal to zero in the first place in the DR scheme calculation. However, this 
naive procedure is technically inconsistent whenever RG running is involved and does not 
work for other calculations involving epsilon scalars, so one should really distinguish between 



the two schemes as a matter of principle. The parameters appearing in the DR effective 
potential obey DR renormalization group equations, which differ from the DR ones with m 2 
set equal to 0. 

The procedure of going from the DR scheme to the DR' scheme is similar at two loops, 
and is described explicitly in section |[ 



4 Two-loop effective potential in the MS scheme 

The two-loop effective potential in the MS scheme for the general theory with masses and 
couplings described by eqs. (|2.6|) - (|2.14 ) can be computed by the methods described in 0. In 



fact, all of the hard work of evaluating the relevant Feynman loop integrals has already been 
accomplished there; no new types of integrals arise. Momentum integrals and vector indices 
each run over 4 — 2e dimensions. For each two-loop diagram, one must include counterterms 
for the various one-loop divergent subdiagrams. The result still includes single and double 
poles in e, which are then simply removed by two-loop counterterms in modified minimal 
subtraction. The final result can be divided into parts corresponding to the various graphs 
of figure Because the VV, VVV, and ggV graphs all involve the same field-dependent 
coupling g abc , it is natural to combine their contributions into a pure gauge piece • 
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For the result, I find: 

y{2) _ y{2) y{2) y(2) y(2) y {2) y (2) y {2) 

v — V SSS + V SS + V FFS + K ^p5 + •'ssv + v vs + ^vvs 

y (2) y (2) y(2 ) 

^ v FFV + FFV + 'gauge 

where, in terms of the masses and couplings as specified in eqs. (|2.6|) -( |04|) , 



(4.1) 



T/(2) _ 


^(X^ffsssimlmlml), 


(4.2) 


y(2) _ 

v ss — 


^fss(mlmj). 


(4.3) 


y(2) _ 
v FFS ~ 


^\y Uk \ 2 fFFs{mj,m 2 j,ml), 


(4.4) 


= 

FFS 


-/'V^M^M^f^simj, m 2 , m\) + c.c, 


(4.5) 


y(2) _ 

v ssv — 


\{g aij ?fssv{mlrnlml) : 


(4.6) 


y(2) _ 

v vs — 


\g aaU fvs{mlml), 


(4.7) 


v vvs — 


\(g ab yfvvs(ml,mlm*), 


(4.8) 


y{2) _ 
v FFV ~ 


^\gf\ 2 fFFv{mj,mlm 2 a ), 


(4.9) 


V— = 

FFV 


lgfgfM n 'M* JJf f mv (ml ml m 2 a ), 


(4.10) 


= 

gauge 


^(g abc ) 2 f g .u g e(m 2 a ,m 2 b ,m 2 c ), 


(4.11) 



in which all indices on the right side are summed over. The loop-integral functions appearing 
here are given by: 



fsss(x,y,z) 
fss(x,y) 
fFFs(x,y,z) 

/ffs( 3; ' Vi z ) 
fssv(x,y } z) 



fvs(x,y) 
fvvs(x,y,z) 



-I(x,y,z), (4.12) 

J(x,y), (4.13) 

J(x,y) - J(x,z) - J(y,z) + (x + y- z)I(x,y,z), (4.14) 

2I(x,y,z), (4.15) 

-{ (—x 2 — y 2 — z 2 + 2xy + 2xz + 2yz)I(x, y, z) + (x — y) 2 I(0, x, y) 
z^ 

+(y - x - z)J(x,z) + (x-y- z)J(y,z) + zJ(x,y)j 

+2(x + y - z/3)J(z), (4.16) 

3J(x,y) + 2xJ(y), (4.17) 

1 {—x 2 — y 2 — z 2 — lftxy + 2xz + 2yz)I(x, y, z) 

Axy L 

+ {x - z) 2 I(0, x, z) + (y- z) 2 I{0, y, z) - z 2 I{0, 0, z) 
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+(z-x-y)J(x,y) + yJ (x, z) + xJ (y, z) } 

+ + \j{v) + 2J ( Z ) -x-y-z, (4.18) 
fFFv(x,y,z) = -{(x 2 + y 2 - 2z 2 - 2xy + xz + yz)I(x,y,z) - (x - y) 2 I(0,x,y) 
+ (x — y — 2z)J(x, z) + (y — x — 2z)J(y, z) + 2zJ(x, y)\ 

+2(-x -y + z/3)J{z) - 2xJ(x) - 2yJ{y) + (x + yf - z 2 , (4.19) 

fFFv(x,y,z) = 6I(x,y,z) + 2(x + y + z)-AJ(x)-4J(y), (4.20) 
f ga u ge (x,y,z) = -£^{(-x 4 -8x 3 y -8x 3 z + 32x 2 yz + 18y 2 z 2 )I(x,y,z) 

+ {y - zf{y 2 + lOyz + z 2 )I{0, y, z) + x 2 {2yz - x 2 )I(0, 0, x) 

+ (x 2 - 9y 2 - 9^ 2 + 9xy + 9xz + 14yz)xJ(y, z) 

+Ax 3 yz + 48xy 2 z 2 + {22y + 22^ - !6x/3)xyzJ(x)} 

+ (x <-> y) + (x <r+ z). (4.21) 

Symmetry factors have been explicitly factored out of eqs. ( |4.2|) -( fll~Tl) , but fermion-loop 
minus signs and other factors associated with the evaluation of the Feynman diagrams 
are contained in the definitions of the functions. The functions obey obvious symme- 
tries: f S ss( x,y,z) and f gajnge (x,y, z) are invariant under interchange of any two of x,y,z, 
while fss{x,y), f FF s{x,y,z), fpp S (x,y,z), fssv{x,y,z), fvvs(x,y,z), f FFV (x,y,z), and 
f-pp V (x,y, z) are each invariant under interchange of x, y. 

The functions involving vector fields contain factors 1/x, 1/y, and 1/z which appear to 
be singular in the massless vector limit. This is due to the appearance in the Landau gauge 
of vector propagators 

1 I r p^ Jjr p 1 ' I p 2 \ 



which give rise to factors 



i \ p 2 + m 2 — it 



1 



(4.22) 



(4.23) 



p 2 {p 2 + m 2 ) m 2 \p 2 p 2 + m 2 J 

in the loop integrals. The massless limits are actually smooth, and arise often in practice. It is 
therefore useful to have explicit expressions for those massless limits that are not immediately 
obvious. Using eqs. ( |2.29j )-( ]2"!n| ), they are found to be: 

fssv(x,y,0) = (x + y) 2 + 3(x + y)I(x,y,0) + 3J(x,y)-2xJ(x)-2yJ(y), (4.24) 
f vvs (x,0,z) = - ^- ^1(0,0, z) + (^-j)mx,z) + ^-J(x,z) 
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+2J(z), 

-3I(0,0, z ) + 7 -J( z )-^, 
0, 



(4.25) 
(4.26) 
(4.27) 



4xy 



{(43x 2 y + A3xy 2 - 7s 3 - 7y 3 )/(0, x, y) 



fgauge(x, 0, 0) 



+ (2y + 7x)x 2 I(0, 0, x) + (2x + 7y)y 2 I(0, 0, y) 

+ (34xy - 7x 2 - 7y 2 ) J(x, y) } + Ax 2 + 4y 2 + —a;?/ 
19 

— ^-[xJ(x) +yJ(y)} + 5[yJ(x) +xJ(y)\, 

59 23 
13x/(0,0, x) — —xJ(x) + —x 2 . 

6 4 



(4.28) 



(4.29) 



All of the functions vanish (by dimensional analysis) whenever all arguments vanish. 

It may also be of interest to see the individual contributions of the Feynman diagrams 



^gauge- These contributions are listed in Appendix A. 

The classic results of Ford, Jack and Jones for the Standard Model are a particularly 
useful special case of those found in this section, with which I have checked agreement. In 
fact, each type of term that can occur in a general model in MS does in fact arise in the Stan- 
dard Model case; no new types of integrals arise, so that the results of eqs. (|4.12| )- (|4.21| ) could 
be inferred from by some forensic combinatorics. Their functions A(x,y,z), B(x,y,z), 
C(x,y), D(x,y,z), E(x,y), T,(x,y), and A(x, y,z) are respectively equal to the functions 
fssv{x,y,z), -f vvs {x,y,z), f vs (x,y), -f FFV (x,y,z), y/xyf TWv (x,y,z), f vv (x,y), and 
—fvvv(x,y,z) given in this section and in Appendix A. [Note that after the published 
errata of ref. 0, a few further minor typographical errors have been recently corrected in 
the eprint archive version.] 

5 Two-loop effective potential in the DR scheme 

In this section, I report the results for the effective potential in the DR scheme. These 
are obtained by keeping all 4 components of each vector field, but performing momentum 
integrations in 4 — 2e dimensions. The difference, compared to the results for MS, can be 
organized in terms of the extra epsilon scalars with multiplicity 2e. Of course, the SSS, SS, 
FFS, and FFS diagrams in fig. || are unaffected by this procedure. Also, the SSV and 
ggV diagrams are unchanged in going from MS to DR, because in those cases all the vector 
indices are contracted with a loop momentum. The VS, FFV and FFV diagrams yield new 



labelled VV, VVV, and ggV in figure even though these can always be combined into 
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contributions, which we can call eS, FFe and FFe, when the vector line in each case is turned 
into an epsilon-scalar line. In the VVS diagram, a non-vanishing additional contribution 
arises only when both vectors are turned into epsilon scalars; call this contribution eeS. In 
the VV diagram, one or both of the vector lines can become an epsilon scalar, yielding 
contributions to be called eV and ee respectively. Finally, in the VVV diagram, any two of 
the vector lines can be turned into epsilon-scalar lines, resulting in a contribution eeV. 

As discussed in section |3|, the couplings of epsilon scalars have exactly the form indicated 
for vectors in eqs. (|2.12|) - (|2.14 ), except that when an epsilon scalar is involved, the rotation 



to the squared- mass eigenstate basis requires rather than N^ v \ This distinction is 
indicated by replacing the vector index a,b,c, . . . by an epsilon-scalar index a,b,c, . . . on 
the couplings. For example [compare to eqs. (|2.15|) - (|2.16|) 1, the epsilon scalar-epsilon scalar- 
vector, epsilon scalar-vector-vector, and fermion-fermion-epsilon scalar couplings are 

9f ef9 N^N^NS\ (5.1) 



9 = 9 J "'"jygjybf -c g , 

g aic = gff'NgNlfNjp, 

gf = g{T b ) K L N?>N$N$. (5.3) 



gf f9 N£N$N%\ (5.2) 



Then the result in the DR scheme can be written 

= VM + V {2) + V® + VB + V® + V® + W + V % (5.4) 

MS £J ' eeb 1 b b e 1 FFe " €e eeV > V / 



where 



Kf = 


\g mi fes(rh 2 a ,m>), 


(5.5) 


v,S = 




(5.6) 


V FFe — 


-\g! J \ 2 fFFe(mlmlm 2 a ), 


(5.7) 


v® = 

FFe 


{gYgfMV'M^f^ml ml ml), 


(5.8) 


V eV — 


\{g hbc ) 2 fev{rnl,ml), 


(5.9) 


v® = 

ee 


\{9^fee{ml,ml), 


(5.10) 


T/(2) _ 
V eeV — 


lig^rfMmlmlml), 


(5.11) 



with the loop functions given by: 

f eS (x,y) = -2xJ(y), (5.12) 
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f F Fe(x, y, Z 
fl?FeK X i Vi Z 

fev(x,y 

fee(x,y 

f eeV (x,y,z 



-2 JO) + x + y + z, 

2xJ(x) + 2yJ(y) - (x + yf + z 2 , 

4 J(x) + 4J{y) -2x-2y- 2z, 

—4xy — 6xJ(y), 

Axy } 



-x 



y 2 + z 2 — 6xy — xz — yz + (6x + 6y — 2z)J(z) 



(5.13) 
(5.14) 
(5.15) 
(5.16) 
(5.17) 
(5.18) 



This completes the result for the two-loop effective potential in the DR scheme. 



6 Two-loop effective potential in the DR scheme 

As explained in the Introduction and in section [3], it is convenient in models of softly broken 
supersymmetry to go to the DR,' scheme. This scheme is defined so that m\ (the difference 
between the squared masses of epsilon scalars and their vector counterparts) does not appear 
in the beta functions of other couplings, or in the effective potential, or in the equations 
relating pole masses to running masses. Starting from the DR results of the previous section, 
I find that this is done at two-loop order by the following parameter redefinition of soft terms 



appearing in eq. (|2.33|) : 



DR 



[m 



DR/i 



16tt 5 



5l2g 2 C{i)ml 



+ 7lt^i Ylk%kl92[C{k) - \ Cil)] 



"T5r' 



A dr' 



C DR + ( 167r 2)2 



ADR 167T 2 



+5{g A C(t) [2S(R) + iC(i) - 6C G ] m 2 }, 

Y^ k ^ k g 2 C(j)m 2 , 

dajm 2 ) 2 1 (9 2 , r Qfo , n u 2 n2 

2 + (16^17 ~ C g]M 

+g 2 d G C G \M\ 2 m 2 + g 2 ^n i3 C(i)m 2 }. 



(6.1) 
(6.2) 

(6.3) 



If there is more than one simple or U(l) group, then each of the correction terms should be 
summed over subgroups, with a different m 2 for each subgroup. The exception is that the 
term 



(6.4) 



involves a double sum over subgroups labeled a, b. Following these redefinitions, the result 
for the full two-loop effective potential turns out to have the same functional form as if one 
naively took the DR result and set m 2 to 0, removing the distinction between and N^ v \ 
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between hatted and unhatted vector squared-mass eigenstate indices on the couplings, and 
between m\ and m\. It is therefore convenient to define functions which combine the effects 
of the 4 — 2e vectors and the epsilon scalars: 

F vs (x,y) = f vs (x,y) + f eS (x,y), (6.5) 

F V vs(x,y,z) = f vvs (x,y,z) + f eeS (x,y,z), (6.6) 

F FFV (x,y,z) = f FFV (x,y,z) + f FFe (x,y,z), (6.7) 

FpFvix^y^z) = fpF V (x,y,z) + fpp e (x,y,z), (6.8) 

^gaugcO^, V, Z) = /gauged, V, z) + feev(x, y, z) + f eeV (z, X, y) + f&v(V, Z, x) 

+fev(x, y) + fev(y, x) + f eV (x, z) + f eV (z, x) + f eV (y, z) + f eV (z, y) 

+fee(x,y)+fee(x,z)+f ee (y,z). (6.9) 

Note that I use F's rather than /'s to distinguish the DR' functions from the corresponding 
MS functions. 



Therefore, the DR two-loop effective potential is given by: 

T/(2) _ 1/(2) , y(2) y{2) y{2) y (2) y(2) y {2) 



+ T/(2) +T/ (2) +T/(2) 
^ v FFV + "ffv + K gauge) 



(2) 



(6.10) 



where now 



y{2) _ 

v sss — 




(6.11) 




^/ S5 (m 2 ,m 2 ), 


(6.12) 


y(2) 

*FFS — 


^l/ J T/FFsK,m 2 ,m 2 ), 


(6.13) 


v^l = 

FFS 


^"V'^M^Mj^/^K, m 2 , m 2 ) + c.c, 


(6.14) 


T/(2) _ 

v ssv — 


^) 2 Wm 2 ,m 2 ,m 2 ), 


(6.15) 


1/(2) _ 


i^F F5 (m 2 ,m 2 ), 


(6.16) 


T/(2) _ 


^(^) 2 ^y5K,m 2 ,m 2 ), 


(6.17) 


T/(2) _ 

v FFV — 


i|^ J | 2 F FW (^,^,m 2 ), 


(6.18) 


V— = 

FFV 


i^/'M^'M^F^K, m 2 , m 2 a ), 


(6.19) 


V® = 

gauge 


^(^ C ) 2 ^gauge(m 2 , m 2 ,m 2 ). 


(6.20) 
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Here f S ss(x,y,z), f ss (x,y), f FFS (x,y,z), fpp S {x,y,z), and fssv(x,y,z) are given by ex- 
actly the same functions as in MS, eqs. ( f4.12|) - ([4.16|) . The new functions are given by: 



Fvs{x,y) 
F V vs(x,y,z) 



3J(x,y) 
1 



(6.21) 



Axy 



x 2 — y 2 — z 2 — lOxy + 2xz + 2yz)I(x, y, z) 



FFV 



[x,y,z) 



+{x - z) 2 I(0, x, z) + (y- z) 2 I(0, y, z) - z 2 I(0, 0, z) 
+(z-x -y)J(x,y) + yJ (x, z) + xJ (y, z) } 

+\j{x) + \j{y), (6.22) 

-{(x 2 + y 2 — 2z 2 — 2xy + xz + yz)I(x, y, z) — (x — y) 2 I(0, x, y) 
z 1 

+ (x — y — 2z)J(x, z) + (y — x — 2z)J(y, z) + 2zJ(x, y)\ 
+2{-x-y + z/3)J{z), 



F FF v\ x i Vi z ) 



6I(x,y,z) 
1 



(6.23) 
(6.24) 



gauge 



x,y,z) 



Axyz 



-x 



lx 6 y - 8x 6 z + 32x 2 yz + 18y 2 z 2 )I(x, y, z) 



+(y - zf{y 2 + IQyz + z 2 )I(0, y, z) + x 2 {2yz - x 2 )I(0, 0, x) 
+ (x 2 - 9y 2 - 9z 2 + 9xy + 9xz + Uyz)xJ(y, z) 
+ {22y + 22z - 4Qx/3)xyzJ(x)} 
+(x <-> y) + (x <-> z). 



(6.25) 



Despite the appearance of x, y, z in the denominators, these functions again all have 
smooth limits for x,y, z —>■ 0. The non-trivial ones are 



F vvs {x, 0, z) 

F vvs (0,0,z) 
F FFV (x,y,0) 



gauge 



[x,y,0) 



-3I(0,0,z) + ~J(z)-~ 
-(x + y) 2 + 2xJ(x) + 2yJ(y), 
^{(43x 2 ?/ + 43xy 2 - 7a; 3 - 7y 3 )I(0, x, y) 
+(2y + 7x)x 2 I{0, 0, x) + (2s + 7y)y 2 I{0, 0, y) 



(6.26) 

(6.27) 
(6.28) 



+ (3Axy - 7x 2 - 7y 2 )J(x, 
25 



2 2 11 
3x + 3y + ~^ X V 



gauge 



(0,0, x) 



(6.29) 
(6.30) 



[xJ(x) + yJ{y)\ + 5[yJ(x) + xJ(y)], 

71 19 

(0, 0, x) xJix) + —x 2 . 

6 4 

Also, it may be of interest to see the contributions from individual graphs to F ga _ nge (x, y, z 
Those are listed in Appendix A. 
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This completes the result for the two-loop effective potential in the DR scheme. These 
are appropriate for use in any softly broken supersymmetric model, including the MSSM. 
Partial results for the MSSM corresponding to the leading contributions proportional to 
asyf, CisUb, and yf, y\ have been given in refs. |!0| and [jT5] . Several illustrative examples 
and consistency checks are done in section Rl 



7 Renormalization group invar iance of the two- loop ef- 
fective potential in softly broken supersymmetry 



In general, the condition for RG invariance of the effective potential is 



l i <t>i 







v = o. 



(7.1) 



(s) 

Here, A/ are all of the running parameters of the model with beta functions [3\j , and 7] are 
the anomalous dimensions of the scalar fields (pi. At one- and two- loop order, this means 



+ 



yaw 9 



tsl) i>i— 

1 d 
, d 



v (o) 



(7.2) 



(7.3) 



In softly broken supersymmetry, I find that the anomalous dimension matrix for scalar 



fields in the Landau gauge and in either DR or DR is 



7> 



7p 

(5.1) j 
i 

(5.2) j 



16n 



2^1 



(S,l)j 



+ 



(167T 



YikiY J - 5lg 2 C{i) 



1 



^Y~Y»Y^ + Y w Y>»ft2C{k)-C<i)] 



+5lg 4 C(i)[S(R) + 2C(i)-~C G }. 



9 



(7.4) 
(7.5) 

(7.6) 



This can be obtained starting from the general results in the MS scheme in ref. |[24|| , and 
then applying the coupling constant redefinitions needed to transform from the MS to the 
DR or DR' scheme p5| . The eigenvalues of this matrix constrained to the subspace of 
the classical scalar background fields give the anomalous dimensions appearing in eqs. ( |7.2j ) 



and (\T3i ). It should be noted that because of gauge-fixing, the Landau gauge scalar field 
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anomalous dimension matrix 7$ relevant for the effective potential is not the same as the 
more widely-known, gauge-invariant, anomalous dimension matrix of the chiral superfields. 



For comparison, the latter is [26 



-A = ^7P + 7^7P; (7.7) 



16tt 2 11 (16tt 2 ) 2 11 
tF = \Y m Y^-28ig 2 C{i) 1 (7.8) 

7P = ~Y mr X nkl Y klr Y mr ^ +Y lH Y^ l g 2 [2C(k)-C(t)} 

+5ig 4 C(i)[2S(R) + 4C(i) - QC G ]. (7.9) 

In order for the effective potential to satisfy eq. ( |7.1| ) in a model with explicit supersym- 
metry breaking, it is necessary to include a running vacuum energy term A, as in eq. ( |2.33| ). 



Now using the results of section |5], one can derive the DR beta function for A up to two loops 
in a general softly-broken supersymmetric model as specified in subsection |2.4j , by looking 
at the ^-independent piece of eqs. (|7.2|) - (|7.3|) - I find 



16tt 2Pa (16tt , 

ffl = (m 2 y i (m 2 y j + 2(m 2 )y k fi k3 + b^b lJ -d G \M\\ (7.11) 
/5| 2) = g 2 d G \M\ A [AS(R) - 8C G ] + 8g 2 \M\ y^-C(z) + 8g 2 {m 2 yJ h ^ kj C{i) 
+4g 2 {m 2 )j{m 2 ))C{i)+4g 2 b lJ b l] C(i)-4g 2 M^b ij C(i)-4g 2 M*^ 
-Y^ k Y l3l \{m 2 )T{m 2 ) l m + (m 2 )™/i mn/ ^ + ^ mn {m 2 ) l n + /i fcm (m 2 )>"' + b km b ml 



[m 2 ) l k + fi kmf i ml - 2Y^ k Y am (m 2 y jf i mn fi nk 



-Y ijk a ijmm b ml - Y ijk a^ l n km b ml . (7.12) 

where d G is the dimension of the adjoint representation of the group. If the gauge group 
contains more than one simple or £7(1) component, then terms involving the gaugino mass 
M or g 2 should be summed over subgroups in eqs. (|7.5|)-(|7T6|), (|7.8|)-([7T9|), and ( [7.11| )-( [TT2j ), 



with the exception noted in eq. ( |6.4| ). Special cases of these general results will be used in 
the next section. 



I have checked explicitly that the DR two-loop effective potential for a general softly- 
broken supersymmetric theory satisfies RG invariance, using the results given above and in 
Appendix B and in refs. [[H], p~0f] . 
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8 Examples and consistency checks 



In this section, I study some examples chosen as consistency checks and useful points of 
reference for the results given above. The examples are all based on supersymmetry with or 



without soft breaking, so the DR scheme is used. One type of consistency check follows from 
the requirement that the two-loop effective potential satisfies RG invariance in conjunction 
with the known two-loop beta functions [ ^T] , |2"2"| , [], ID], and the scalar anomalous dimensions 



and P\ found in the previous section. (Since the result for /3a is itself a consequence of 
the calculation, that part is not really an independent check.) The derivatives of the loop 
functions are listed in Appendix B. Another type of check relies on the fact that the effective 
potential for a supersymmetric theory in a supersymmetric vacuum must vanish. These 
consistency checks rely on non-trivial cancellations between different two-loop functions, 
which are made manifest by writing them in terms of the basis functions I(x,y, z), J(x,y), 
and J(x), using eqs. (P^)-(P^) and (|OTD-(|O0D. 



8.1 The Wess-Zumino Model 



Consider the Wess-Zumino model |27| with a single chiral supermultiplet <£> containing a 
Weyl fermion ip and a complex scalar <j> + (R + il) / v2, where 4> is the classical background, 
and R, I are real scalar quantum fluctuations. The superpotential is given by 

W = ^ + H^ (8.1) 
2 o 

where /i and y are mass and coupling parameters, taken to be real without loss of generality. 
The fields R, I, ip are mass eigenstates, with 

m \ = fi 2 + 3yfi(f) + 3y 2 (f) 2 /2, (8.2) 
mj = fi 2 + yii<j ) + y 2 <j ) 2 /2, (8.3) 
rrty = fi + y<p. (8.4) 

The tree-level scalar potential is 

V<® = /iV + + y 2 4 /4, (8.5) 

and the one-loop contribution to the effective potential is given in terms of the function h(x) 
in eq. ( |3.13| ) by 

V {1) = h(m 2 R ) + h(mj) ~2h(ml). (8.6) 
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The non-zero scalar quartic and cubic couplings are: 



^RRRR _ y 

yiRii 
ytRR 
ytu 



,iiii 



and the Yukawa interactions are 



V 

y 



i>i>R 



y 2 l% 

3y(M + #)/A 
y(/i + #)/V2, 



//A 



;s.7) 
;s.8) 
;s.9) 

5.10) 



^ = iy/V2. 



(8.11) 
(8.12) 



It follows that the contributions to the two-loop effective potential are: 

,2 



Vc 



(2) 



V 



sss 



(H + y4>) 3f SS s(m R , m R , m R ) + fsss(m R , m z , m I 



(2) 



ss 

(2) 
FFS 

(2) 
FFS 



y_ 

16 

y 2 



3fss(m R , m R ) + Sfssirrij, m 7 ) + 2f ss (m R , m I 



r I 2 2 2\ , t I 2 2 2 

jFFsym^ m^m R )+ fpFsKm^, m I 



y 2 2 



^, uu^ iio r j -r JFFSyn^, 

r / 2 2 2 \ r ( 2 2 2\ 

JFFS\ m i)i m i)i m R) JFFS\ m -tpi m ^i m Ij 



Now one may check RG invariance of the effective potential. At one-loop order, 
from eq. ( |S.6|) that 

Q^ (1) = -M 2 -^ 3 -^ 4 /4. 

The one-loop scalar anomalous dimension and beta functions are 



(8.13) 
(8.14) 
(8.15) 
(8.16) 
one finds 

(8.17) 



y 2 /2, 
y 2 v, 

3y 3 /2. 



(8.18) 
(8.19) 
(8.20) 



Therefore, from eq. 



d 



2t/V 2 2 + 5y 3 /i0 3 /2 + 3t/ 4 4 /4, 
-y W - 3y 3 ^ 3 /2 - y^/2, 



(8.21) 
(8.22) 
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where A/ runs over y, /z, so that eq. ([7.2|) is indeed satisfied. At two-loop order, one finds 
from eqs. (|8~T3|)-(|8TT6|) and that 



(i) 
A/ 



9 



,(5,1). 



5 



yd) 



0Aj ' c 

From the two-loop RG scalar anomalous dimension and beta functions: 



(8.23) 



7?< 2) 



4 

3y 5 /2, 



one also finds: 



9 



2/V 2 + 3y 5 W>72 + y 6 072. 



(8.24) 
(8.25) 
(8.26) 

(8.27) 
(8.28) 



The results of eqs. ( |8.23|) , ( |8.27|) and ( |3.28| ) combine to verify eq. (|?T3|). 

In the special case of (ft = 0, supersymmetry is not broken, and the effective potential 
should vanish. At one-loop order, eq. ( |S.6|) then vanishes trivially. At two-loop order, 

,2 

<.2<)) 



V{2) = \ [^fsSsiH 2 ,^ 2 ,^ 2 ) + fss{H 2 ,V 2 ) + /W 5 (/i 2 ,// 2 ,/i 2 ) 



which equals by virtue of eqs. (|4.12|) - ([4.14j ). 



8.2 Supersymmetric QED in supersymmetric vacua 

Let us now consider a supersymmetric U(l) gauge theory with coupling constant g and a 
pair of chiral superfields with charges ±1. 

First take the case that the chiral superfields do not have a mass term before symmetry 
breaking, and the two scalar fields have the same classical background value (ft. Then the 
gauge symmetry is broken, but supersymmetry remains unbroken since (ft parameterizes a 
flat direction. The vector boson, two Weyl fermions, and a real scalar field each obtain a 
mass 



x = Ag 



2 J,2 



(8.30) 



Together with a massless (in Landau gauge) real scalar Nambu-Goldstone boson, these form 
a massive vector supermultiplet. In addition, there are two massless real scalars and one 
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massless Weyl fermion forming a chiral supermultiplet. The DR one-loop effective potential 
vanishes because of these mass degeneracies. The two-loop effective potential contributions 



in the DR scheme are: 



v sss 


= 9 2 x ^-fsss(0,0,x) + fsss(0,x,x) , 
12 


(8.31) 


V FFS 


= 9 2 [fFFs(0, x, 0) + f FFS (0, ^> x) + 2f FFS (x, X, 0)] , 


(8.32) 


v ssv 


q 2 

= j[fssv{0,0,x)+f SS v{0,x,x)), 


(8.33) 


v vs 


= y Fvs ( x,x " >} 


(8.34) 


v vvs 


= g 2 xF vvs (x,x,0), 


(8.35) 


v (2) 

v FFV 


= g 2 F FFV (Q,x,x), 


(8.36) 



with the other contributions vanishing. One can now check by plugging in the results of 
section |6] that the sum of eqs. flB.31 )-( 8.36| ) yields 0, as required for a supersymmetric vacuum. 



This constitutes a non-trivial identity involving cancellations between different two-loop 
functions which become apparent after writing them in terms of the functions I(x,y,z), 
J(x, y) and J(x). 

Another check which relies on a different set of cancellations is obtained if we take = 
in the above model, but now include a superpotential mass term \x. In that case, the vector 
gauge boson and the gaugino are massless, and the real scalar fields and the chiral fermions 
all have squared mass /i 2 . Then one obtains for the contributions to the two- loop effective 



potential in the DR scheme: 

V® = 9 2 fss(^^ 2 ), (8.37) 

V F % = VW0,/rV), (8.38) 

Vssv = 9 2 fssv(^,^,0), (8.39) 

V F 2) FV = g 2 F FFV (fi 2 ^ 2 ,0), (8.40) 

Vg y = -sVi^^/AO), (8.41) 

with all other contributions vanishing. Again one finds from the results of section || that the 
sum of eqs. ( |8.37|) - (|8.41|) yields 0, as required for a supersymmetric vacuum. 
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8.3 Supersymmetric SU(N C ) gauge theory with one flavor in su- 
persymmetric vacua 

A richer set of checks is found in non-abelian supersymmetric models. As an example, 
consider supersymmetric SU(N C ) gauge theory with one flavor of chiral superfields Qi and 
Q in the fundamental and anti-fundamental representations, respectively. Here i = 1, . . . , N c 
is a color index. Consider evaluation of the effective potential for the classical background: 

(Qi) = (Q 1 ) = <W- (8.42) 

These VEVs break the gauge symmetry according to SU(N C ) — > SU(N C — 1), but </> param- 
eterizes a flat direction and supersymmetry is unbroken. Therefore the effective potential 
must vanish at each order in perturbation theory for any value of <fi. My aim is to show this 
explicitly. 

The particle content for non-zero (f> consists of 2N C — 1 massive vector supermultiplets with 
their associated massless (in Landau gauge) real scalar Nambu-Goldstone modes, N 2 — 2N C 
massless vector multiplets associated with the unbroken gauge symmetry, and one massless 
singlet chiral supermultiplet. The non-zero squared-mass eigenvalues are 



x 



9 , 
2(iV c 



,2 i2 



(8.43) 
(8.44) 



and the multiplicities of the mass eigenstates are shown in Table |I[ Because of the mass 



Table 1: Multiplicities of mass eigenstates in the model of section |0 



particle type 


m 2 = 


2 

m = 


- X 


m 2 = y 


real scalars 


2N C + 1 


2N C - 


- 2 


1 


Weyl fermions 


N 2 C -2N c + l 


4iV c - 


-4 


2 


vectors 


N 2 - 2N C 


2N C - 


- 2 


1 



degeneracies indicated in Table [3], the one-loop contribution to the effective potential vanishes 
as required. 



At two-loop order, I find the contributions in the DR scheme to be 



K 



(2) 



-4/2 



SSS 



^lf sss (0, 0, x) + {Nc m p 2 fsss(0, 0, y) + ^-^fsss(0, x, x) 
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v ss 

V FFS 



V-, 



(2) 



FFS 



(2) 



ssv 



(2) 



vs 



(2) 

vvs 



W 2) 



FFV 



(2) 



y(2) 



gauge 



, (iVc-l) 2 f m , (iV c - 2) 2 (iV c - 1) 
+ ^2 /sss(0, y, y) + ^2 /sss(0, x, y) 







9 



'2JV? - 3iV c - 1 



5.45) 
(8.46) 



{/ FFS (0, x, 0) + / FFS (0, x, x)} + 3(iVc 2 ^ fFFsix, x, 0) 



TV — 1 
N 2 r - N c + 2 



v fFFs{y,y,Q) 

c 



N P -1 



9 



{fFFsfa 2/, 0) + f F Fs(x, y, x)} + iVc 2 x fFFs(x, x, y) 

2 {f~FFs( X ' X i ^ ) _ f~FF~s( X > X i ^)} 



Af — 1 1 

/W (0, 0, x) + -fssv(0, 0, y) + 



4Ar JsSv{0,y,V) 



+ 



4 

AL -1 



-f 



J5sy(0, x, x) H fssv{x, x, 0) + —fssv{x, x, y) 

fssv(x, y, x) 



4 

' — - — F y5 (x,x) + — - — F y <j(x,y) + -^-F vs {y,x) 
N c -1 



4iV, 



-Fvs{v,y) 



9 



A A? 



N C (N C - 2) 



7V„- 1 



(iV c - 2) 



F yy5 (0, x, 0) + -^_F yy5 (x, x, 0) 
(iV c - I) 2 



S—^F vvs {x, y, 0) + 1 C N2 ' F vvs (y, y, 0) 



2 A" 2 -3N C -1 



N c -l 



FpFV 

(0,x,x) + N c (N c -2)F FFV (x,x,0) 



^-^F FFV (0, y, y) + ^/V ~ FpFV ( x > x ' ^) + ~\~ F ' 



FFv(x,y,Xj 



-9 



Aj3 



N C (N C - 2)F ¥Fv (x, x, 0) + Fppyix, x, y) 



+2(N c -l)F FFv (x,y,x 
~N C (N C - 2) 



-fgauge (0, X, x) + i^g aU g e (x, X, ?/) 



5.47) 



(8.48) 



(8.49) 



5.50) 



(8.51) 



3.52) 



(8.53) 



(8.54) 



After some algebra, using eqs. ( |4.12|) -( flT6"|) and ( |6.21|) -( |O0|) , one finds that the sum of these 
contributions indeed vanishes, as required by unbroken supersymmetry in the flat direction 
parameterized by 4>. 
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8.4 Softly-broken supersymmetric QED 

Consider the case of supersymmetric QED with a coupling g and two chiral superfields with 
charges ±1, as in subsection fi.2\ However, now we introduce supersymmetry-breaking effects 
in the form of a gaugino mass M, and non-holomorphic soft supersymmetry-breaking scalar 
squared masses m 2 + and m 2 for the scalar fields of charge +1, —1 respectively. Instead 
of equal VEVs, the scalar fields of charge +1, —1 are taken to have classical background 
values <ft, respectively. Then the four real scalar mass eigenstates obtain squared masses 
Xi, Xi, x 2 , x 3 where 



Xx 



2 2/2 

m_- g , 



2 i 2/2 

x 2 = m + + g <p , 
x 3 = m 2 + 3g 2 4> 2 , 



(8.55) 
(8.56) 
(8.57) 



and the three fermion mass eigenstates obtain squared masses 0,yi,y2, with 



yi = [M 2 + 4g 2 <j) 2 - M^/M 2 + 8g 2 <f) 2 ]/2, 
y 2 = [M 2 + Ag 2 (j) 2 + MJM 2 + 8g 2 (j) 2 ]/2, 



(8.58) 
(8.59) 



while the vector boson obtains a mass 



z = 2g 



2 j. 2 



(8.60) 



Because supersymmetry is explicitly broken, RG invariance requires that a vacuum-energy A 
is included among the soft supersymmetry breaking terms. The tree-level potential is then: 

,2 



V® =A + 7740 2 + |-0 4 . 



From eq. ( |3.13| ), the DR one-loop effective potential contribution is: 

V {1) = 2h( Xl ) + h(x 2 ) + h(x 3 ) - 2h( Vl ) - 2h(y 2 ) + 3h(z). 



(8.61) 



(8.62) 



In that scheme, by following the procedures described in sections |2.1| and || I find the 
following contributions to the two-loop effective potential: 



(2) 



sss 

(2) 



ss 



r 1 3 

g 4 4> 2 [fsss{xi,xi,x 3 ) + -fsssfa, x 2l x 3 ) + -fsss(x 3 ,x 3l x 3 ) 

r 1 1 3 

9 2 [fss(xi,xx) - -fss(xi,x 2 ) - ^fssfai, x 3 ) + -fss{x 2l x 2 ) 

1 3 

+ ^fsS (X 2 ,X 3 ) + -fss(x3,X 3 ) 



(8.63) 



(8.64) 
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(2) 



FFS 



2(?/i + z)f FFS (0, yi, xi) + 2(2/2 + z)f FFS (0, 2/2, xi) 



M 2 + 8g 2 (j) 2 

+2z{f FFS (yi, yi, x 2 ) + fFFs(Vl, Vu x z) + fFFsfa, 2/2, x 2 ) + fpFsfa, 2/2, 2:3)} 



FFS 



^yi{fTFs(y^y^ x ^) - fFFsivuvu^)} 



+M 2 {f FFS (y u 2/2, x 2 ) + f FF s(yu 2/2, x 3 )} 

M 2 + 8g 2 (f> 2 
+22/2{/ffs(?/2, 2/2, a*) - f FF s(V2, 2/2, ^2)} 

+M 2 {f ¥Fs (y l , 2/2, x 2 ) - fpp 3 (yi, 2/2, x 3 )} 



(8.65) 



V, 



(2) 



ssv 



(2) 



9_ 
2 

„2 



vs 

(2) 

vvs 



fssv(xi, xt, z) + fssv(x 2 , x 3 , z) 



F vs (z,x 1 ) + ^F vs (z,x 2 ) + ^F vs (z,x 3 ) 



g 2 zF vvs (z,z,x 3 ), 



v FFV 



FFV 

gauge 



(M + 8g 2 <p 2 )F FFV (V } 0, z) + y 2 F FFV (y u 2/i, z) 



2(M 2 + 8c/20 2 

+yiF FFV (y 2 , 2/2, 2) + 2zF FFV (y h 2/2, « 

22/ 6 4 
M 2 + 8c/ 2 02 

0. 



F ffv(v^V^ z) + Fpp V (y 2 , y 2 ,z)- 2F FFv (y 1 , y 2 , z) 



(8.66) 

(8.67) 

(8.68) 
(8.69) 

(8.70) 
(8.71) 
(8.72) 



We can now test the RG invariance of the effective potential. The one-loop scalar anoma- 



lous dimension and beta functions in the DR scheme are: 



Pm 



7f" 

a 



(i) 

M 
(1) 



-g , 
2g\ 

4g 2 M, 

-8g 2 M 2 + 2g 2 (m 2 + 
-8g 2 M 2 + 2g 2 (m 2 _ 



m 



m 



+ h 



(ml) 2 + (m 2 ) 2 - M A . 



(8.73) 
(8.74) 
(8.75) 
(8.76) 
(8.77) 
(8.78) 



From eq. Q8.62 ) one therefore finds that 

M 4 - (m 2 + ) 2 - (m 2 ) 2 + 8g 2 (P 2 M 2 + 2g 2 <p 2 m 2 _ 



- W™+ - AgY: 



(8.79) 



and, from eqs. (g7g)-(g7g 



(!) 9 y(0) 



M 4 + (ml)' 2 + (m 2 _)' 2 - 8g 2 <f) 2 M 2 - 2g 2 4> 2 m 2 



d\i 
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+2g 2 (f) 2 m 2 , + 2g 



4 ±4 







so that 



2<? 2 2 m 2 + 2<? 



4/4 



so that eq. (|7.2|) is satisfied. At two loop order, one has 



lf 2) 
Pm 



(2) 







1 m_ 
Iff 



V, 
8g 5 , 

32 g 4 M, 

96g 4 M 2 + 16g 4 m 2 + , 
96g 4 M 2 + \Qg A m 2 _ , 
V(m 2 ) 2 + Ag 2 (m 2 _ f + 8# 2 M 4 



V/? (2) — O) 



8<? 2 M 4 + Ag 2 {mif + % 2 (m 2 _) 2 + 96^ 4 2 M 2 



+16g 4 (j) 2 ml + 8g 



6^4 



One also finds from eqs. ( ^.63| )- (|8.72| ) and the results of section ||: 

y{l) 



lg A (j> 2 m 2 + - 8#V- 







9A/ 



JO) 
51) 



^2) 
53) 
54) 
55) 
?6) 
57) 



(8.89) 



g 2 M 4 — 4g 2 (m 2 + ) 2 — Ag 2 {m 2 _) 2 



-9Qg*M 2 (f) 2 - 8# 4 <fm: 



4 j.2,2 



.90) 



Together, eqs. (|8~88D -(|8^0D verify eq. Q. 



9 Outlook 



In this paper, I have presented the results for the two-loop effective potential of a general 



renormalizable field theory in the Landau gauge, in each of the MS, DR, and DR renormal- 
ization schemes. These results should be useful in connecting specific models of electroweak 
symmetry breaking to future data in a precise way. 

It is not unlikely that the correct model for physics near the TeV scale is based on some 
version of softly-broken supersymmetry, either the MSSM or some moderate extension of it. 
Previous calculations of the effective potential in the MSSM have used the one-loop result 
and partial two-loop approximations with leading corrections proportional to asyf and 
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yf [29|-|33|. However, there is still some RG scale-dependence in these results, compared to 
estimates of our eventual ability to measure properties of the Higgs sector at future colliders. 
Use of the full two-loop DR' effective potential should further reduce the scale dependence. 



RG improvement methods [I7|-[20|, P^|-||38|| should enable an accurate determination of the 
vacuum of the MSSM and its extensions. I plan to report on the application of the results 
of the present paper to the MSSM soon 



Appendix A: Individual diagram contributions to the 
functions / gauge and F gauge 

The three Feynman diagrams labelled VV, VVV, and ggV in figure |2| all involve the field- 
dependent coupling g abc , and combine to yield V^ ge . In the MS scheme, the individual 
diagram contributions to the function f gavi ge{x,y, z) are given in an obvious notation by 

f ga ,u g e(x,y,z) = fvvv(x,y,z) + f vv (x,y) + fvv(x,z) + fvv(y,z) 

+ f 99 v(x) + fggv(y) + fggv(z), (A.l) 

where 

f vvv (x,y,z) = -^—{(-x 4 -8x 3 y-8x 3 z + 32x 2 yz + 18y 2 z 2 )I(x,y,z) 
Axyz L 

+(y - z) 2 (y 2 + lOyz + z 2 )I(0, y, z) - x 4 J(0, 0, x) 

+ (x 2 — 9y 2 — 9z 2 + 9xy + 9xz — 13yz)xJ(y, z) 

„ s 129 2 n /20x y z\ 
+Axyz + — xy z - I — + - + - 1 xyzJ[x) j 

+ (x <-» y) + (x <-> z), (A. 2) 

fvv(x,y) = 2 ^J(x,y) + ^J(y) + ^J(x) + 6 -^-, (A.3) 
fggv(x) = |/(0,0,x) + |j( i r). (A.4) 

Similarly, in the DR' scheme, 

F gauge (x,y,z) = F vvv (x,y,z) + F vv (x,y) + F vv (x,z) + F vv (y,z) 

+f 99 v(x) + f 99 v(y) + fggv(z), (A.5) 



where 



F vvv (x,y,z) = fvvv(x,y,z) + f eeV {x,y,z) + f eeV (z,x,y) + f eeV {y,z,x), (A.6) 
Fvv(x,y) = f vv (x,y) + fev(x,y) + fev(y,x) + f ee (x,y). (A.7) 
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with fggvix) given as before. Explicitly, 

F vvv (x,y,z) = -^^{(-x 4, -8x 3 y -8x 3 z + 32x 2 yz + 18y 2 z 2 )I(x,y,z) 

+{y - z)V + Wyz + ^ 2 K(o, y, z) - x 4 /(o, o, x) 

+ (x 2 — 9y 2 — 9z 2 + 9xy + 9xz — 13yz)xJ(y, z) 

xy 2 z 2 . AAx My 47z. T . ,1 
+ (-— + -f + — )xyzJ(x)\ 

+ (x <-> y) + (x <-> z), (A. 8) 

Fvv(x,y) = 2 ^J(x,y)-fj(y)-fj(x)-^. (A.9) 



The results for vanishing arguments are easily obtained from eqs. (|2.29|) - (|2.31 



Appendix B: Renormalizat ion-group-scale derivatives 

It is often useful to have expressions for the derivatives of the two-loop effective potential 
functions with respect to the renormalization scale Q, for example to check RG invariance. 
The derivative of the one-loop effective potential function h(x) defined in eq. ( 3.1 3|) is 

Q^h(x) = -x 2 /2. (B.l) 

For checking the RG-invariance of the effective potential, it is convenient to write the deriva- 
tives of two-loop functions with respect to Q in terms of the derivative of the one-loop 
function with respect to its squared-mass argument: 

h'(x) = ^(Ex - 1). (B.2) 

•Li 

The derivatives of the two-loop functions can all be found from those of the basis functions: 

Q-^Ax) = -2x, (B.3) 

Q-^J(x,y) = -Ayh'(x)-Axh'(y) J (B.4) 

Q-^I(x,y,z) = A[h'(x) + h'(y) + h'(z)]-2(x + y + z). (B.5) 
For the derivatives of the MS two-loop functions, one finds: 

Q-^fsss(x,y,z) = -A[h > (x) + h'(y) + h'(z)}+2(x + y + z) J (B.6) 

Q-^fss(x,y) = -4yh'{x)-Axh'(y), (B.7) 
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Q-^f F Fs(x,y,z) = 4xh'(x)+4yh'(y) + (8x + 8y-4z)h'(z) 

-2x 2 - 2y 2 + 2z 2 - Axy, (B.8) 

Qj^f-FFS&y,*) = S[h'(x) + h'(y) + h'(z)]-4(x + y + z), (B.9) 

Q-^fssv(x,y,z) = 12xh'(x) + 12yh'(y) + (12x + 12y-Az)h'(z) 

-2x 2 - 2y 2 - 12xy - 6xz - 6yz + ^-z 2 , (B.10) 

Q-^fvs(x,y) = -12yh'(x)-12xh'(y)-Axy, (B.ll) 

Q-^fvvs(x,y,z) = -9h'(x)-9ti(y)-12ti(z) + ^ + ^-z, (B.12) 

Q-^f F Fv(x,y,z) = (12x + 12y-8z)h'(z) + 6xz + 6yz + — : (B.13) 

Q-QQfFFv(x,y,z) = 24[h'(x) + h'(y) + h'(z)}-4x-4y-12z, (B.14) 

Q-Q^f ga ngc(x,y,z) = 9[(x + y)h'(z) + (x + z)h'(y) + (y + z)h'(x)\ 

+52[xh'(x) + yh\y) + 
19 

-^-(^ + I/ 2 + - 63(xy + yz + a;*). (B.15) 
For the functions used with epsilon scalars in the DR scheme, one has: 

Q^fes(x,y) = ±xy, (B.16) 

Q-^Us(x,y,z) = 4z, (B.17) 

Q^f F Fe(x,y,z) = -Ax 2 -Ay 2 , (B.18) 

Q-^}fFFe(x,y,z) = -8x-8y, (B.19) 

Q^fev(x,y) = I2xy, (B.20) 

Q-^fee(x,y) = 0, (B.21) 

Q^feev(x,y,z) = -12xz-12yz + 4z 2 . (B.22) 

Finally, the functions used in the DR,' scheme (besides those found in MS) satisfy: 

Q-^F vs (x,y) = -12yh'(x)-12xh'(y), (B.23) 

Q-^F vvs (x,y,z) = -9h'(x)-9ti(y)-12h'(z) + — + f + 3z, (B.24) 
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Q—F FFV (x,y,z) = (I2x + Uy-Sz)h'(z) -Ax 2 - Ay 2 + 6xz + Qyz + — , (B.25) 
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